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We introduce a local formalism, in terms of eigenstates of number operators, having 
well defined point symmetry, to solve the Hubbard model at weak coupling on a TV x N 
square lattice (for even N). The key concept is that of W — states, that are the 
many-body eigenstates of the kinetic energy with vanishing Hubbard repulsion. At 
half filling, the wave function demonstrates an antiferromagnetic order, a lattice step 
translation being equivalent to a spin flip. Further, we state a general theorem which 
allows to find all the W = pairs (two-body W = singlet states). We show that, in 
special cases, this assigns the ground state symmetries at least in the weak coupling 
regime. The N = 4 case is discussed in detail. To study the doped half filled system, 
we enhance the Group Theory analysis of the 4x4 Hubbard model introducing an 
Optimal Group which explains all the degeneracies in the one-body and many-body 
spectra. We use the Optimal Group to predict the possible ground state symmetries of 
the 4x4 doped antiferromagnet by means of our general theorem and the results are 
in agreement with exact diagonalization data. Then we create W = electron pairs 
over the antiferromagnetic state. We show analitycally that the effective interaction 
between the electrons of the pairs is attractive and forms bound states. Computing 
the corresponding binding energy we are able to definitely predict the exact ground 
state symmetry. 



I. INTRODUCTION 

The (repulsive) 2D Hubbard hamiltonian is one of the popular models of the high-Tc7 cuprates [1], 
as many people believe that it contains at least some of the relevant ingredients of the mechanism of 
superconductivity. While other ingredients may well be missing for the full explanation of supercon- 
ductivity and the rich phase diagram of these materials, there are now strong evidences of pairing in 
this model, and although pairing is not synonimous to superconductivity it can hardly be supposed 
to be extraneous to it. The evidence for pairing comes from various independent methods, including 
cluster diagonalizations [2] [3] [4], fluctuation exchange (FLEX) [5] diagrammatic approach, which is 
based on a conserving approximation, and renormalization group techniques [6] [7] . The approach we 
are proposing is based on an analytic canonical transformation, reminiscent of the original Cooper the- 
ory. However, our understanding of the many-body ground state is only partial, and the symmetries 
of the bound pairs have not yet been fully explored. The Hubbard model must be more thoroughly 
understood before we can solve more realistic ones. 

In the strong coupling limit the double occupation of the same site is energetically suppressed and 
the model at half filling is equivalent to the Heisenberg model with an antiferromagnetic exchange 
interaction [8]. A popular approach takes care of the strong repulsion between two opposite spin 
fermions by a Gutzwiller [9] projection, i.e. by throwing out of the Hilbert space the double occupation 
states. However, truncating the Hilbert space in this way costs kinetic energy, so at finite U the system 
must allow double occupation, also in the ground state, as one can see from the eigenvectors of cluster 
calculations. At weak coupling, on the other hand, it makes sense to speak about particles in filled 
shells, which behave much as core electrons in atomic physics, and particles in partially filled, or 
valence, shells. Remarkably, particles in partially filled shells can totally avoid double occupation at 
no cost in energy; they do so, forming W=0 states, that are defined as many-particle eigenstates 
of the kinetic energy with no double occupation. Below, using a new formalism, we show how 
W = states arise by symmetry. We stress that since W = states emerge from symmetry alone, 
they remain W = for any coupling strength and are an adequate starting point for a realistic theory. 
This is the reason why weak coupling expansions often provide good approximations at intermediate 
coupling, as observed by several Authors [10] [11] [12]. 
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A. Canonical Transformation Approach to the Pairing Mechanism 

Let us consider the Hubbard model with hamiltonian 

H = H + W = tJ2J2 C r* C r'° + 12 Ufl rinrl, U > 0, (1) 

cr {r, r ') r 

on a square lattice of N x N sites with periodic boundary conditions and even N. Here a =|,| 
is the spin and r, r' the spatial degrees of freedom of the hole creation and annihilation operators 

and c respectively. The sum on (r, r') is over the pairs of nearest neighbors sites and h ra is the 
number operator on the site r of spin a. The point symmetry is C\ Vl the Group of a square [13]; 
besides, H is invariant under the commutative Group of Translations T and hence the Space Group 
[14] G = T <g> C± v ; (g> means the semidirect product. We represent sites by r = (i x ,i y ) and wave 
vectors by k — (k x , k y ) — jf(i x , iy), with i x , i y = 0, . . . , N — 1. In terms of the Fourier expanded 
fermion operators Cfe CT = ^2 r e %kr c ra , we have H = J2k e (^) c ka c k<r wl tli e (k) = 2t[cosk x + cosfc y ]. 
Then the one-body plane wave state cj.jO) = \ka) is an eigenstate of Ho- 

To study the behaviour of two holes added to the system in its ground state, we introduce the 
W = pairs, a special case of W = states. Using degenerate eigenstates of the kinetic energy Ho, 
the non-abelian point symmetry Group C iv of the Hubbard hamiltonian (1) allows the existence of 
two-body singlet states with no double occupancy. They are obtained by a configuration interaction 
mechanism and may be of special interest since the particles of a W — pair do not feel any direct 
repulsion. Hence a new pairing scenario arises: since the two extra particles cannot interact directly, 
by definition of W — pair, their effective interaction comes out from virtual electron-hole excitations 
exchanges with the Fermi sea and in priciple can be attractive. In the following we show how to get 
the effective interaction between two holes added to the system. 

Many configurations contribute to the interacting (n+2)-body ground state \^o{n+2)) and we need 
a complete set S to expand it exactly; as long as it is complete, however, we can design S as we please. 
We can take the non- interacting n-body Fermi sphere \<f>(n)) as our vacuum and build the complete 
set in terms of excitations over the vacuum. In the subspace with vanishing spin z component, the 
simplest states that enter the configuration mixing are those obtained from 1 3>(n)) by creating two 
extra holes over it (we call them the m states); the set 

4 lT 4 2 J $ ( n ))' <k 1 ),e(k 2 )>e F (2) 

where e_p is the Fermi energy, is a basis for this part of the Hilbert space. Similarly, along with the 
pair m states, we introduce the 4-body a states, obtained from |$(n)) by creating 2 holes and 1 
electron-hole (e-h) pair; a basis looks like 

c LT c L4 Cfc ^ c L<rl $ W>' e(k 1 ),e(k 2 ),e(k 4 ) > e F > e(fc 3 ), (T=T,| (3) 

Then S includes the 6-body (3 states having 2 holes and 2 e-h pairs, and so on until states with two 
holes and n e-h pairs if we are below half filling or 2N 2 — n e-h pairs if we are above. We are using 
Greek indices for the configurations containing the electron-hole pairs, which here are playing largely 
the same role as phonons in the Cooper theory. By means of the complet set S we now expand the 
interacting ground state 

|*o(« + 2)) = J2 a mH + J2 ba ^ + Jl^ll 3 ) + ■- ( 4 ) 

m a (3 

and set up the Schrodinger equation 

H\* (n + 2)) =E\Ko(n + 2)). (5) 

We stress that Eq.(4) is configuration interaction, not a perturbative expansion. 

When the number n of holes in the N x N system is such that |$(n)) is a single non-degenerate 
determinant (the Fermi surface is totally filled), we can easily and unambiguously define and calcu- 
late the effective interaction between the two extra holes since the expansion (4) for the interacting 
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ground state is unique: this is done by a canonical transformation [15], [16], [17] from the many-body 
Hamiltonian of Eq.(l). We consider the effects of the operators on the terms of \^o(n + 2)). We write: 

H \m) = E m \m), (6) 

and since W can create or destroy up to 2 e-h pairs, its action on an to state yields 

W\m) = YWm>,m\m') + £>>W a , m + ^|/3)^, ro . (7) 

m' a f3 

The matrix elements of the two-body interaction W between determinants which differ by two spin- 
orbitals reduces to a bielectronic integral. Hence one can show that if to', to represent different W = 
pairs added to the vacuum, W m >^ m vanishes. Moreover, 

W\a) = ^\m)W m , a + Xy>W a ,, Q + Y)®Wfi, a + h)W ya , (8) 

ma' j3 7 

where scattering between 4-body states is allowed by the second term, and so on. In this way we 
obtain an algebric system for the coefficients of the configuration interaction (4). However to test 
the instability of the Fermi liquid towards pairing it is sufficient to study the amplitudes a m of the 
to states. In the weak coupling limit this can be done truncating the expansion (4) to the a states 
because, as we have shown [15], the inclusion of the (3, 7, . . . states produces a renormalization of the 
matrix elements of higher order in W, leaving the structure of the equations unaltered. Choosing the 
a states such that 

(H + W) aa , = 5 aa ,E' a (9) 

the algebric system reduces to 

(E m - E)a m + ^a m -W m , m / + ^2b a W„ ha = (10) 

m' a 



(E' a -E)b a + Y,a m ,W a , m , =0. (11) 

m' 

Solving for b a and substituting in the first equation we exactly decouple the 4-body states as well, 
ending up with an equation for the dressed pair \a) = J2 m a ml m )- The effective Schrodinger equation 
for the pair reads 

{H + W + W eff [E])\a)=E\a) (12) 

where 

(W ef f) mm > = - 2^ E , _ E ■ (13) 

a a 

Hence Eq.(12) is a self-consistent equation and we have calculated the effective interaction W e ff 
between W — pairs analytically [16], [17]; it can be attractive depending on n. Also we have found 
that the results compare well with exact diagonalization results, when available. Basically the same 
mechanism works for small clusters with open boundary conditions [3]. 

We want to stress that the truncated expansion of | x f'o(^ + 2)) in Eq.(4) doesn't give a good approx- 
imation of the interacting ground state wave function but only of its weak coupling a m amplitudes. 
Similarly in the BCS model, from the Cooper equation (obtained by truncation of \^o{ n + 2)) to the 
to states) we can estimate only the pair coefficients of the ground state and not its full structure. 
Nevertheless this is enough to study bound states formation; indeed the energy gap of the pair in the 
Cooper theory and in the many-body BCS theory are equals. 
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B. Truncated Configuration Interaction: Description of the Pairing Mechanism in Small 

Clusters 

Although the above canonical transformation can be performed without any truncation of the 
configuration interaction expansion, the description simplifies if we take into account a states only. 
As observed above this is justified in the weak coupling regime and is enough to study the pairing 
problem. Let us consider for example the pairing mechanism in small clusters, like C11O4 and Cu^Oiq. 
Such clusters with 4 holes have ground states of 1 B 2 symmetry in G± v for < U < t. The diagnosis 
that hole pairing between the holes of a W = pair occurs in these ground states results from the 
following steps. 

a) The lowest eigenstate of the one-body energy spectrum has A\ symmetry and the interacting 
ground state with two holes is a d v totalsymmetric singlet. The state |<I>) = 0)4^0^ JO) is the two- 
body U — non-degenerate ground state, where here and in the following steps c\ a is the creation 
operator for the hole eigenstate of I symmetry and spin a. b) The first excited level of the one-body 
energy spectrum is degenerate and the corresponding eigenstates transform as the vector components 
(x, y). c) To have a ground state with four holes of 1 B 2 symmetry in the limit of vanishing interaction 
we must create two holes over |$) in the singlet state \ms) = ^75 ( c lt c yi. + c yT c xO^- ^ Taking \ms) 
as our unique m state we apply the above canonical transformation and we find an effective attractive 
interaction summing over all the virtual 4-body (3 holes-1 electron) intermediate states. 

As another description to this weak-coupling case, we may say we are doing degenerate second-order 
perturbation theory; in this example, the zero-th approximation eigenstates are just two, \m\) = 
cL.cLj$) and \rri2) = c^*cLj<l>), an d the problem is reduced to a 2 x 2 matrix diagonalization. The 
corresponding second-order interaction is illustrated in Figure I: it takes the state |mi) into \m 2 ) and 
conversely. Hence it is actually a spin-flip interaction. 




Figure I. The e-h exchange diagram for the Two-Hole Amplitude. For W = pairs, the direct interaction 
vanishes and this diagram produces an effective interaction, splitting singlet and triplet pairs. 

It turns out [3] that the effective interaction is attractive in the singlet \^ib 2 ) = \ms) = ^(\ m i) + 
\m 2 )) and equal in size, but repulsive, in the triplet |^1/3 B2 ) = |m<r) = -t=(|mi) — |m,2})- That is 

definitely a pairing situation for the singlet, e) The two-body state ^= ( c lt c li + C j/T c ^l^^ is a PF = 
pair (no direct interaction occurs between the holes); hence we can say that the W = pair, when 
dressed by virtual e-h excitations, gives up a bound pair. 

These are the principles that here we wish to extend to the half filled 4x4 Hubbard model of 
Eq.(I). 



C. Plan of the paper 

Below, we consider the ground state of the 4x4 cluster with 14, 15 and 16 holes, and demonstrate 
that the 14 hole case must be interpreted in terms of pairing of two electrons added to the half-filled 
antifcrromagnet. However the non- interacting ground state is degenerate and this means that we must 
extend the above formalism. 

The physical interpretation emerges from the fact that we are able to solve the problem analytically 
at weak coupling. When the analytic results are compared to the numerical ones [2] we can conclude 
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that i) we are able to predict from Group Theory the good quantum numbers and degeneracies of the 
ground states involved at various occupancies ii) we are able to show that the effective interaction 
between the extra electrons is attractive in the singlet (ground state) and repulsive for the triplet. 

These analytic results have been allowed by a deeper symmetry analysis of the 4x4 cluster than 
had been possible previously. The antifcrromagnetic ground state has also been explored analytically 
by a new approach. In this way, the electron pair creation is also accoplished analytically with full 
control of the symmetry. 

Due to the special role played by W — pairs, in the next Section we state a general theorem to 
obtain all the possible W — pairs. The theorem requires the knowledge of a (previously unknown) 
symmetry Group big enough to explain the one-body degeneracies (Optimal Group). As we shall see 
the Space Group G does not work and in Section III we determine the Optimal Group for the 4x4 
square lattice. In Section IV we determine the exact weak coupling ground state wave function at 
half filling, that is unique as granted by the Lieb's theorem [18]. This will be done using a new local 
formalism that enables us to write down all the W = states of the partially filled shell. In Section 
V, with the help of the theorem of Section II and with the Optimal Group in hands, we single out all 
W = pairs formed from degenerate orbitals of the shell e(k) = 0. Hence in Section VI we use the 
corresponding local basis to write down W — two and four-body states. Then to study the pairing 
problem at half filling we extend the above canonical transformation and finally in Section VII we 
shall reduce the pairing problem to 2 x 2 matrices and compare with the literature numerical data. 
Finally, we underline the implications of the present results in Section VIII. 

II. W=0 PAIRS AND SYMMETRIES OF THE DOPED GROUND STATES: A USEFUL 

THEOREM 

In general, one has a set iS of degenerate eigenstates of Ho which is partially filled in the U = 
limit. To first-order, we ignore the particles in the filled shells and find the exact ground state(s) of 
W in the truncated Hilbert space 7i describing those in the partially filled shell. One simple case 
occurs when only two holes are in the partially filled shell, because the ground states are W = pairs 
(filled shells are understood). We recall that they are two body singlets that are eigenstates of Ho and 
belong to the kernel of W. There is no double occupation of any site in such a pair. This means that 
the two particles of a W = pair do not interact directly, but only by means of virtual electron-hole 
excitations. 

As we shall see, another special case occurs when the system at half filling is doped by removing 
two particles (they can only be removed as W = pairs). 

In this Section we want to show how W = pairs arise, due to the symmetry of the system. 
Previously [15] [16] we have shown that W = pairs with zero total momentum are a consequence of 
the Ci v point symmetry. Projecting the determinantal state 



on the irreducible representations A 2 , B\ and B 2 of C\ v one obtains W — pairs. 

Here we want to point out a more powerful and elegant criterion to get all the W — pairs, 
including those of nonvanishing total momentum. We can do that in terms of the Optimal Group Q of 
the Hamiltonian, that we define as a symmetry Group which is big enough to justify the degeneracy 
of the single particle energy levels. By definition, every one-body eigenstate of H can be classified 
as belonging to one of the irrcps of Q. We may say that an irrep r\ is represented in the one-body 
spectrum of H if at least one of the one-body levels belongs to 77. Let £ be the set of the irreps of Q 
which are represented in the one-body spectrum of H . Let \ip) be a two-body eigenstate of the kinetic 
energy Ho with spin S z — and the projection operator on the irrep r\. We wish to prove the 

W=0 Theorem: 



In other terms, any nonvanishing projection of \ip) on an irrep not contained in E, is an eigenstate of 
Hq with no double occupancy. The singlet component of this state is a W — pair. Conversely, any 
pair belonging to an irrep represented in the spectrum must have positive W expectation value. 



|d(fc)>=cU-fcilO> 



(14) 



rj £ £ <^WP M \ip) = 



(15) 
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In the case of the 4x4 model (see Section V) the pairs belonging to irreps of the Optimal Group 
have well defined parities under particle exchange, and can be classified as singlet or triplet pairs. 

Proof: Let us consider a two body state of opposite spins trasforming as the z-th component of the 
irrep r\ of Q: 

l^ } > = £^(^2)4^4^0). (i6) 

Then we have 

n r tn rl \4 v) ) = ^" ) (r,r)4 T 4jO) = (r,r)\r T,r |). (17) 
We define P^ as the projection operator on the i-th component of the irrep r\. Since 

£ ^ ^ OH, r |> = £ 4 V) (r, r)\r T, r |) (18) 

r r 

if 

ff'V T,r|) =0 Vr (19) 
then ^^(r, r) = Vr. It is worth to note that Eq.(19) is true if and only if 

where \ra) — c\. a \Q). 

It is always possible to write \ra) as 

i^) = EE c i 7?) wbS) ( 2 o) 

r)££ i 

where \^pfl) is the one-body eigenstate of Hq with spin a transforming as the i-th component of the 
irrep r\. From (20) it follows directly that if r( does not belong to £ 

pM\ra) = 

and so pW)\ r |, r |) = 0. □ 

This theorem restricts the possible ground state symmetries in some special case. Let U c (n) be the 
minimum crossover value of U, that is the ground state with n holes \^o{n)) has well defined symmetry 
770 for < U < U c (n). Let |<&(n)) be non-degenerate (closed shells case). Then if we add two extra 
holes, the new ground state for U — is a W — pair over |$(n)), and its symmetry is 770 ■ i]w=o, 
where r)w=o is the symmetry of the added pair. Turning on the interaction the symmetry cannot 
change ii U < U c (n + 2)! This restriction is posed by Group theory alone: which of the symmetries 
that remain allowed is actually realised in the ground state depends on the dynamics. 

The complete characterization of the symmetry of W — pairs requires the knowledge of the 
Optimal Group Q. A partial use of the theorem is possible if one docs not know Q but knows a 
subgroup, like the Space Group G. It is then still granted that any pair belonging to an irrep of G 
not represented in the spectrum has the W — property. On the other hand, accidental degeneracies 
occur with a subgroup of the Optimal Group, and by mixing degenerate pairs belonging to irreps 
represented in the spectrum one can find W — pairs also there. This is clearly illustrated by the 
example repported in the next Section. 



III. THE OPTIMAL GROUP FOR THE 4x4 HUBBARD MODEL 

The half filled shell of the 4x4 Hubbard model has degeneracy 6, and the Space Group G does 
not have irreps with dimensions bigger than 4 [14]. An additional symmetry is needed to justify the 
sixfold degenerate eigenvalue e(k) — 0, and we have found it [19]. Let us represent the 4x4 lattice as 
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1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 



Periodic boundary conditions arc assumed and for example, the nearest neighbours of 1 are 2, 5, 4 
and 13. Rotating the plaquettes 1,2,5,6 and 11,12,15,16 clockwise and the other two counterclockwise 
by 90 degrees we obtain the effect of the "dynamical" symmetry, that we call d: 



5 


1 


4 


8 


6 


2 


3 


7 


10 


14 


15 


11 


9 


13 


16 


12 



This transformation preserves nearest neighbours (and so, each order of neighbours) but is not an 
isometry, and for example the distance between 1 and 3 changes. Thus, this symmetry operation d is 
a new, dynamical symmetry. Including d and closing the multiplication table we obtain the Optimal 
Group Q with 384 elements in 20 classes (like G) as shown in Table I. 



Ci 


c 2 


c 3 




c 5 


c 6 


c 7 


c s 


c 9 


C10 


I 


t 22 


da' [2] 


Ox 


c 2 


a' 


C 2 d 


C 2 t 22 d 


C 2f r'[l] 


C%to 2 d 


C u 


C12 


C13 


C14 


Cl5 


Ci6 


Cl7 


Ci8 


C19 


C 2 Q 


C$t 2 od 




cwi] 


c 4 


C 2 to\d 


C 2 t 12 d 


C 2 a x [l}d 


C 4 [l]d 


C 2 [l]d 


C 4 [l] 



Table I. //ere, we report one operation for each of the 20 classes d; the others can be obtained by conjugation. 
The operations are: the identity I, the translation t mn of m steps along x and n along y axis; the other opera- 
tions C2, C4, a, a' are those of the Group of the square and are referenced to the centre; however,C2[i], Ca[i\, a[i], 
and a'[i] are centered on site i. 

The complete Character Table of Q is shown as Table II. 
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Q 


Ci 


c 2 


c 3 


C 4 


c 5 


c 6 


c 7 


c 8 


C 9 


C10 


Cn 


C12 


Cl3 


Cl4 


£-15 


Cx6 


C 17 


Cis 


£-19 


C20 


A x 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


Ai 


1 


1 


1 


1 


1 


1 


-1 


-1 


1 


-1 


-1 


-1 


-1 


1 


1 


1 


1 


1 


-1 


-1 


B 2 


1 


1 


-1 


-1 


1 


1 


1 


1 


1 


-1 


-1 


1 


-1 


-1 


1 


1 


-1 


-1 


1 


-1 


Bo 


1 


1 


-1 


-1 


1 


1 


-1 


-1 


1 


1 


1 


-1 


1 


-1 


1 


1 


-1 


-1 


-1 


1 


1 


2 


2 


_2 


_2 


2 


2 








2 














_2 


-1 


-1 


1 


1 










2 


2 


2 


2 


2 


2 








2 














2 


-1 


-1 


-1 


-1 








Si 


3 


3 


3 


3 


3 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 














1 


1 


S9 


3 


3 


3 


3 


3 


-1 


1 


1 


-1 


1 


1 


1 


1 


-1 














-1 


-1 




3 


3 


_3 


_3 


3 


-1 


-1 


-1 


-1 


1 


1 


-1 


1 


1 














1 


-1 


s 4 


3 


3 


-3 


-3 


3 


-1 


1 


1 


-1 


-1 


-1 


1 


-1 


1 














-1 


1 


Ai 


4 


-4 


-2 


2 








-2 


2 





-2 


2 











-1 


1 


-1 


1 








A 2 


4 


-4 


-2 


2 








2 


-2 





2 


-2 











-1 


1 


-1 


1 








A 3 


4 


-4 


2 


-2 








-2 


2 





2 


-2 











-1 


1 


1 


-1 








A 4 


4 


-4 


2 


-2 








2 


-2 





-2 


2 











-1 


1 


1 


-1 








fii 


6 


6 








-2 


-2 


-2 


-2 


2 








2 


























^2 


6 


6 








-2 


-2 


2 


2 


2 








-2 


























3 


6 


6 








-2 


2 








-2 


-2 


-2 





2 

























6 


6 








-2 


2 








-2 


2 


2 





-2 























ni 


8 


-8 


-4 


4 
































1 


-1 


1 


-1 








n 2 


8 


-8 


4 


-4 
































1 


-1 


-1 


1 









Table II. Character Table of the Optimal Group Q of the 4x4 model. 



As the notation suggests, the irreps A\ and A\ both reduce to A\, in d Vl while B 2 and B 2 both 
reduce to B 2 - Table III shows how the irreps of Q split in C& v . 



Q 


Civ 


A x 


A x 


Ax 


A x 


B 2 


B 2 


B 2 


B 2 


Ti 


2B 2 


r 2 


2Ax 


Si 


Ax + 2B X 


s 2 


Ax + 2B X 


S3 


2A 2 + B 2 


s 4 


2A 2 + B 2 


Ai 


Ax+Bx+E 


A 2 


Ax+Bx+E 


A 3 


A 2 + B 2 + E 


A 4 


A 2 + B 2 + E 


Qx 


A 2 + Bx + 2E 


n 2 


A 2 + Bx + 2E 


n 3 


Ax+B 2 + 2E 


n t 


Ax+B 2 + 2E 


iii 


2Ax + 2Bx + 2E 


n 2 


2A 2 + 2B 2 + 2E 



Table III. Reduction of the irreps of Optimal Group Q of the 4x4 model rn the point Group. 
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We call G the Optimal Group because it enables us to explain the degeneracy of the one-particle 
energy spectrum; in other terms, no accidental degeneracy of orbitals occours using Q. In the Table 
IV below we report the one-body eigenvalues for t = — 1, the degeneracy and the symmetry of 
each eigenvector. Below, we shall find that Q is also adequate to classify the many-body ground 
states. For pairs, the W = Theorem ensures that no double occupancy is possible in the irreps 
A x , B 2 , Ti, T 2 , Si, £ 2 , S 3 , E 4 , A 2 , A 3 , (l u Q 2 , n 3 , H and U 2 . 



Energy 


Irrcp of 


Degeneracy 


4 


S 2 


1 


2 


A 4 


4 







6 


-2 


Ai 


4 


-4 


Ax 


1 



Table IV. One-body spectrum for t = — 1. 

IV. LOCAL FORMALISM AT HALF FILLING 

As observed in the previous Sections, W = pairs predict the possible ground state symmetries of 
systems which differ from closed shells by a pair. Let us now consider how the above analysis extends 
to the doped half filled system. Let Shf denote the set (or shell) of the k wave vectors such that 
e(k) = 0. At half filling (N 2 holes) for U — the ground state has the Shf shell half occupied, while 
all |fc) orbitals such that e(fc) < are filled. The k vectors of Shf lie on the square having vertices 
(±7r, 0) and (0, ±7r); one readily realizes that the dimension of the set Shf, is \Shf\ = 2A — 2. 

For N = 4, the 6 wave vectors are k\ = (tt, 0), & 2 = (0,7r), k% = (tt/2, ir/2), k± = (it/2, — ir/2), k$ = 
(— 7r/2, -7t/2) and fc 6 = (-tt/2, tt/2). 

(0,tt) 




Figure I. The Brillouin Zone; the dashed square marks the condition of vanishing kinetic energy. Shf con- 
tains 6 states (arrows) belonging to one- dimensional irreps of T. Moreover, the arrows are mixed by op- 
erations of C^v, and (tv, 0) and (0, n) are the basis of a two-dimensional irrep of G, while the wavevectors 
k — (±7r/2, ±7r/2) mix among themselves and yield a four- dimensional irrep. The degeneracy of these two 
irreps is accidental for G, but is explained by the Optimal Group Q. 
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Since H commutes with the total spin operators, 




s~ = (s+y 



(21) 



r r 



2 



at half filling every ground state of H is represented in the S z = subspace. Thus, H has 



degenerate unperturbed ground state configurations with S z — 0. We will show how this degeneracy 
is removed by the Coulomb interaction W already in first-order perturbation theory. Actually most of 
the degeneracy is removed in first-order, and with the help of Lieb's theorem [18] we shall be able to 
single out the true, unique first order ground state of H. In Appendix A we show that the structure 
of the first-order wave functions is gained by diagonalizing W in the truncated Hilbert space Ti 
spanned by the states of 3 holes of each spin in Shf- In other terms, one solves a 6-particle 
problem in the truncated Hilbert space Ti. and then, understanding the particles in the filled shells, 
obtains the first-order ground state eigenfunction of H in the full 16-particle problem. We underline 
that the matrix of Hq in H is null, since by construction Tt is contained in the kernel of Ho. 

The operator ^ r h r ^h r i has eigenvalues 0,1,2,... and so the lowest eigenvalue of W is zero (in 
other terms, W is positive semi-definite). The unique ground state of the Hubbard Hamiltonian for 
U = + at half filling will turn out to be a W — singlet state of 6 holes in Shf (filled shells being 
understood). We shall obtain the W — states £ H. It is clear that, although the U = case is trivial, 
at U = + we are still facing a bona fide many-body problem, that we are solving exactly [20]. In the 
present Section we define a basis of local orbitals for the 4x4 Hubbard model with periodic boundary 
conditions; this basis is crucial for making the problem tractable, both at half filling and for the doped 
case. Using the local basis, the many body wave function of the antiferromagnetic ground state can 
be projected out as the singlet component of a single determinant, which is amazingly simple for an 
interacting system; the effective interaction between the doped holes also emerges analytically. The 
treatment for the half-filled case has already been generalized [21] to the N x N Hubbard model; 
below, we present the much simpler solution of the 4x4 case, which is sufficient for our present 
purposes. 

Since W depends on the occupation number operators h r , it is intuitive that its properties in H 
are best discussed by a suitable one-body basis of Shf such that at least one of these operators is 
diagonal. In addition, a convenient basis should exploit the large G symmetry of the system. If Sh / 
were a complete set (iV 2 = 16 states), one would trivially go from plane waves to atomic orbitals by 
a Fourier transformation; instead, we must define the local counterparts of plane-wave states using 
only the 2 N — 2 = 6 states that belong to Shf ■ 

For each site r we diagonalize the number operator n r (for the moment we omit the spin index); 
it is a trivial matter to verify that (n r )ij = (ki\n r \kj) = -L e A ki ~ kj ) r has eigenvalues 3/8 and five 
times 0. This degeneracy suggests that we should diagonalize other operators in order to label the 
h r eigenvactors, and indeed, since h r is compatible with the operations of the point symmetry group 
C^v we also diagonalize the Dirac characters of this Group. The set of Dirac characters defines the 
irreducible representation (irrep); thus we write the one-body basis states {| where a comprises 

the h r eigenvalue and an C\ v irrep label. It is easy to verify that for r — the eigenvector with 
nonzero eigenvalue is just the totally symmetric superposition of all the {\h}} € Shf- Translating by 
r, plane wave states pick up a phase factor: — ► e lkr \k). Hence the h r eigenvector of occupation 



3/8 is 




(22) 



and we set up our local basis at r by 



6 




(23) 



6 




(24) 
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where we introduce the orthogonal matrix 



O = 





- 1 


1 


1 


1 


1 


1 ■ 




V2 


V2 


-1 


-1 


-1 

V2 


-1 

V2 


1 








fi 


V 2 


[J 

V 2 


[J 

V 2 




V3 


-V3 



























-V3 
















V3 





-Vs. 



(25) 



It is clear from Eqs.(23)(25) that \<j)^) has well defined occupation n and symmetry for point group 
operations centered at site r; namely, it has n = 3/8 for i = 1 and n — otherwise; it belongs to Ai 
for i = 1, 2, to B2(xy) for i = 3, to B\(x 2 — y 2 ) for i = 4 and to i? for i = 5, 6. The local bases of 
different sites r and r/ are connected by the unitary transformation 



ICO = 



(26) 



and using the orthonormality of the |fc) states we obtain the elements of the symmetric translation 
matrix 



->(r',r) 



) ^ ^ Oj,m.Oi,rn& 



■k m .(r—r') 



(27) 



m— 1 



The translation matrix knows all the G symmetry of the system, and must be very special. Using 
such a basis set for the half filled shell the antifcrromagnetic order of the ground state comes out in a 
clear and transparent manner. It is clear that [T^ r ,r )] 4 = 1. The translation by one step towards the 
right is accomplished by 



/ji(right) 








-1 

V3 



The matrix that makes one step upwards is 



-i(lip) 



1 

2 



VS V3 

—i —i 

Ve Ve 

i —i 



V2 





i 

VS 


V6 


V% 










V3 


— i 

V6 


—t 
V2 










s is 





















i 

VI 


i 
V3 


—i 
VS 











pi 

V 3 


—i 
V& 


i 

Ve 














i 
V2 


75 


1 

V3 


V 3 














i 

VS 


— i 

V§ 


i 
^/2 











V3 


76 


75 












V2 







(28) 



(29) 



The reason why this choice of the basis set is clever is now apparent. The local basis at any site r 

splits into the subsets S a = {\<f>i'^), |02^)> I03^)}i an< ^ — {|04^)j |05^)i l^^}}] a shift by a lattice 
step sends members of S a into linear combinations of the menbers of Sb, and conversely. 
Consider the 6-body determinantal eigenstate of H a 



M (r) [l,2,3]) c 



(30) 



the notation implies that |e^ r )[z, j, k]) a denotes a 6-body determinant with one body per local state 
and k with spin a, the complement with spin — cr; local states are ordered in the natural way 
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1, .. .6. In this state there is partial occupation of site r with spin a, but no double occupation. 
Introducing the primitive translation of the lattice e x = (1,0) and i y = (0, 1), it turns out that a shift 
by a lattice step r — ► r' = r ± e ; with I = x,y, produces the transformation 

|d< r >[l,2,3]>„ — |^') [1,2,3]), = -|rfW[4,5,6]) CT , (31) 

that is, a lattice step is equivalent to a spin flip (antiferromagnetic property). Since the spin- flipped 
state is also free of double occupation, |d( r ) [1, 2, 3]) CT is a TT = 6-body eigenstate of iJ. A ground 
state which is a single determinant is a quite remarkable property of an interacting model like this, 
and this property holds at half filling, not in general. To be sure, \S r ^ [1, 2, 3]) CT is a mixture of pure 
spin components with 5 = 0,1, 2, 3. However, W is positive semi-definite and this implies that 

all the pure spin components must possess the W — property as well. In particular, the singlet 
I^af ) is a TT = eigenstate and is the true ground state of the Hubbard model at half filling which 
is predicted by Lieb's theorem (filled shell are understood). Explicitly, the antiferromagnetic ground 
state wave function reads 

\$ s A f) = i(3,6)i(2,5)i(l,4)|rfW[l,2,3]) CT , (32) 

where A is the antisimmetrizer, such that for example A(l, 4)|c£( r ) [1, 2, 3]) CT = \d^ [1, 2, 3}) a — 
|d«[4,2,3]) CT . One can easly verify that 1$^°) i s independent by the r and a label of [1, 2, 3]) CT , 
modulo phase factors. From the general analysis of Ref. [21] we obtain that this singlet has A\ sym- 
metry with respect to the center of an arbitrary plaquette of the square lattice and vanishing total 
momentum as in the strong coupling limit [22] . It is worth noticing that the open-shell part of the 
antiferromagnetic ground state (not considering the occupied inner shells) is a 6-body W = singlet 
state. Correlation effects enable no fewer than 6 particles to completely avoid double occupation in 
such a small system. This is also a consequence of Lieb's theorem. If all the 6 body are taken with 
parallel spin, double occupation is trivially avoided; however, Lieb's theorem enforces a singlet ground 
state, so a singlet W = state must exist. 

The 4x4 case at hand can be thorougly explored on the computer, since the size of Ti at half filling 
is 400. We have used Mathematica to diagonalize H + <;S 2 , where a small £ is a numerical device 
to keep the different spin components of the ground state separated. In this way, we observed the 
fourfold degenerate, W = ground state which £ separates into its singlet, triplet, quinted and septet 
components, as expected. At £ = the separation grows like U 2 . The antiferromagnetic property of 
the wave functions was also easily and nicely borne out by the numerical results. 



V. W = PAIRS AND QUARTETS IN THE PLANE- WAVE REPRESENTATION. 

In this Section we use the antiferromagnetic ground state I'&^p ) to predict the possible symmetries 
of the doped half filled system. With 12 holes, in the U — > limit, there are two in Shf, the first-order 
ground states correspond to W = pairs. The symmetry of these W = states can be determined a 
priori from the Q irreps of Table IV. Apart from the filled shells, two holes go to the ^4 level (Table 
IV). From the character Table II one can derive that 

Sll = Ai + B 2 + 4 + Ti + T 2 + £ 2 + S 3 + Si! + 2 + ft 3 ; (33) 

since the first 3 entries are present in Table IV, the W = Theorem ensures that Ti, T2, S2, S3, fii, S^2 
and f2 3 pairs have no double occupation. It turns out that the spin and orbital symmetries are 
entangled, i.e. some of these pairs are triplet and the rest singlet. We can see that by projecting 
the determinantal state cj^cL |0) with k,p e Shf on the irreps not contained in the spectrum. One 
obtaines singlet W = pairs for Ti, T 2 , £2, ^li they read 

i^j ri ) = {^ c lAi + i T 4a) - ^(4 3 t4 3 i - clAi + 4 5T 4 5 , - 4 6T 4a)}|o> 
i^ ri > = 3(4 8 t4.i - 4 4T 4 6i + 4 sT 4 3 i - 4 6T 4a)i°> (34) 

\^^ 2 ) = {ts(4it4ij. + 4 2 t4 2 i) + 2vi(4 3 t4 5 4 + 4 4 t4 6 j. + 4 5 t4 3 j. + 4 6 t4 4 i)}I°) 
iv^ r2 ) - 3(4st4 8 i + 4 4 t4 4 i + 4 sT 4 5 i + 4 6T 4 6 i)l°) (35) 
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i^ S2 ) = ^(4 lT 4a-4 2 T c L)i°) 

= fclAi + c h c U + 4 5T 4 6i + 4 6T 4a)i°) ( 36 ) 

1^3 2 ) = 2( C fe 3 T C fe6l + C ki] C k 5 l + C k 5 ] c kil + C k 6 l C k 3 i)\°) 



and finally 



and the two sixfold sets 



and 



iV'i 2 ) = \( c k 1 -\ c k 3 i ~ 4 3 t4i4 + 4 2 t4 5 i - 4 5 t4 2 |)|o) 

1^2 2 ) = 2( C feiT C fe4l _ C k l \ C k 1 \. ~ C k 2 ] C k 6 l + C k^ c k 2 [)\®) 
1^3 2 ) = 2( C feiT C fe5l ~ C fcsT C *:il + C k 2 ] C k 3 i ~ % T C k 2 1 ) I ^) 
1^4 2 ) = 2( C fciT C feei _ C k 6 ] C k 1 l ~ C k 2 t C k 4 l + % T % J. ) I ^) 

i^ n2 ) = ^(4 3 T4 5i -4 5 T4 3 t)io) 
l^> = ^(4 4T 4 6 t-4 6T 4 4 t)|o> 



iV'i 3 ) = i( c fc 1 t4 3 i _ c fc 3 T c Aii - 4 2 t4 5 i + 4 5 t4 2 |)|o) 

1^2 3 ) = 2( c k 1 T; c k i l ~ c fc 4 T c feil + c fc 2 T c *=6l ~ c fc 6 T c fe 2 4-'l^ 
1^3 3 ) = 2( c feiT c fcsi _ c fc 5 T c feil _ c fc 2 T%l + c fc 3 T c fc 2 l 

i^? 3 > = Mt4 k , - cjLtcL + 4,t4, - 4t4,)|o) 



(37) 



iv^ 1 ) = i(4 lT 4 3 t + 4314a + 4 2T 4 5i + 4 5T 4 2i )lo) 

1^2 *) = 2( C feiT C fe4| + C fc 4 T C *:i4 _ C fc 2 T C fcei _ C *=6 T % J. ) I ^) 
1^3 X ) = 2( C feiT C fe 5 | + C /csT C *:i4 + C fc 2 T%l + C /£3T C '=2i''l^ 
1^4 *) = 2( C feiT%i + C feeT C *:il ~~ C k 2 ] C kil ~ C /c 4 T % J. ) I ^) 

i^ ni ) = 7f(43T43i-4 5 T4 5i )i°) 
i^ 1 ) = 73(4 4T 4 4i -4 6T 4 6i )i )- 

Other irreps yield W = triplet pairs. They are the three times degenerate irrep 

l^ S3 ) = 7f(4 lT 4 2i "4 2 T4a)|0> 

i^ 3 > = §(4 sT 4a - 4 4T 4 3i + 4 5T 4 6i - 4 6T 4a)i°> ( 38 ) 

1^3 3 ) = 2( C fe 3 T%l ~ C k^ C ksl + C fe 5 T C fe44 ~ C k^ c k 3 i>\^) 



(39) 



(40) 



The above a priori argument hardly applies to the symmetries of W — quartets, because f2| contains 
almost every symmetry and we do not know any W — 0-like Theorem for quartets. However, we can 
still build the projection operators by Mathematica; we can project the 225 quartets on the irreps of Q 
and carry on the analysis in an efficient, if not elegant, way. We found that the singlet W = quartets 
are 13 as many as the singlet W ~ pairs, and belong to the same irreps T\, I^, £2, Therefore, 
these are the possible symmetries of the first-order ground states with 14 holes. Exact diagonalization 
results [2] show that for U/t < 3 and 16 — 2 = 14 holes the ground state is sixfold degenerate, with 
a doublet of states with momentum (tt 7 it) and a quartet with momentum (±7r/2, ±ir/2). In view of 
Eqs.(37), the computed ground state corresponds to an Qi electron pair over the half filled system. 
For U/t > 3 and the same number of holes a level crossing takes place: the ground state is threefold 
degenerate and contains a state with momentum (0, 0) and a doublet with momentum (it, 0) and 
(0, 7r). In view of Eqs.(36), the computed ground state must be assigned to a X2 electron pair over 
the half filled system. In both cases, the symmetry of the ground state corresponds to a W — pair. 
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VI. LOCAL PICTURE OF THE W = PAIRS AND QUARTETS 



The antiferromagnetic property of the local basis of any site r readily prompts W — pairs; 
alternatively, we can transform the pairs and quartets of well-defined symmetry (Section V) using 
the local basis. Both methods are useful. The vacuum at half filling is a W — state with 6 holes. 
The Quartets are then obtained from the antiferromagnetic ground state at half filling by removing a 
W = pair. 



A. W=0 Pairs 

Chosen a site r, let \<pa^) be a normalized linear combination of the states \4>^}, |</ > 2*' ) )> 1^3*^) °f S a 
and <j)^ be a normalized linear combination of \4>^), \4>£^), |</>6^) °f Then, 

|d£i) = \4$*H) ( 41 ) 

is a two-body state free of double occupation on every site and ^^ /3T - is a VF = pair. 

Since for each spin one has 3 degrees of freedom, one can build 9 independent pairs in this way; they 
are bases for the fii and £2 W — pairs of the previous Section. The only alternative method for 
obtaining W — pairs is that of forming \4>a^} and \<t$P) as linear combinations of states of the same 

subset (both from S a or both from Sb). This can be accomplished in such a way that |</>^) never 
appears for both spins; actually, the bases of T\ and are obtained in this way. If we use S a for both 
spin directions then after a lattice step in any direction the pair is formed exclusively with states of 
Sb and its occupation vanishes; such pairs live on a sublatticc. 

We rewrite the pairs and quartets of Section V using the local basis of any site r. By Eq.(24), 

MM>= E e- i{ki+ki)r OniO m M { 3<t>i\)- (42) 

m,n— 1 

For instance, the operator that annihilates the 6-th component of the Oi W = pair in Eq.(37) 
becomes 

1 1 1 1 

^ 6 ni = ~^( c iT c 6| + C6TC14) + 2( C 6T C 3| + c 3 TC6|) + ^=(c 2T c 6i + c 6T c 2i ), (43) 

(r)^ (r) . ... . (r) 



where c\J and c l(T are hole creation and annihilation operators in the local states 4>\J ; in Eq. (43) the 
site is not specified since whatever it is \I/ 6 1 does not change. The local representation of symmetry- 
adapted pairs is of interest because different irreps are well characterized by their local behaviour. 



B. W=0 Quartets 



With 14 holes, in the t7 — > limit, there are four in Shf, the first-order ground states correspond 
to W = quartets. By removing in all possible ways two holes of opposite spin from the 6-body 
W = determinant \d^ [1, 2, 3])<j of equation (30) one produces nine 4-body determinants. They are 
free of double occupation on site r because 10^) cannot appear for both spins; this property holds 
on all sites because of the special form of the translation matrices (28) and (29). It follows that these 
are W — states and in first-order perturbation theory they belong to the ground state multiplet. 
It is clear that these 9 first-order ground states are in one-to-one correspondence to the pairs (41), 
and for instance | ^24 V5 j ) corresponds to the pair \(f)^ (f)^}; since |$^°) transforms as the 
totalsymmetric one-dimensional irrep A\ under the C^ v operations referred to the center of a plaquette, 
the symmetries of the 9 quartets and the 9 pairs are also the same , namely, they are fii and £2 of Q . 
The total momcntun labels are also the same. There is actually a complete correspondence between 
W = pairs and quartets; the quartets are also 13 and those that were not obtained above belong 
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to Ti and T 2 . These cannot arise in the same way because one can show that they are not obtained 
by removing two holes from \S r ^ [1, 2, 3]) CT . This means that Ti and T 2 are not to be interpreted as 
pairing states. 



VII. PAIRING MECHANISM 

We consider the ground state of the 4x4 model with 14 holes; aside from the 10 holes in the inner 
Ai and Ai shells (see Table IV) the outer shell contains 4 holes in a W = quartet. We are in 
position to show that the principles of Section I B produce a diagnosis of pairing between two electrons 
added to the antiferromagnetic 16-holcs ground state (half filling). 

We recall from Section V that, by comparing with exact diagonalization results [2], the ground state 
is assigned to fli, Eq.(37), at weak coupling and to X 2 , Eq.(36), at a stronger coupling. In analogy 
with the points a)-c) of Section IB, we must preliminarily verify that symmetry does not forbid 
obtaining these symmetries by creating W = electron pairs, Eqs. (34-40), from the antiferromagnetic 
state of Eq.(32). This is the same as annihilating hole pairs. Since the state l^^p ) of Eq.(32) is a 
total symmetric singlet with vanishing total momentum, the labels of the quartets will be the same 
of the annihilated hole pairs. This operation can be done by hand, or with the help of Mathematica, 
and the answer to the preliminary question is adfirmative for Cl\ and X 2 , but not for all pairs. We 
obtain 24 W — quartets out of the 28 states in Eqs. (34-40), since the annihilation of Ti and T 2 
W = pairs gives identically zero. 

We need to modify the canonical transformation to deal with the 4x4 cluster with periodic boundary 
conditions near half filling; the previous form is not adequate because the Shf shell is only partially 
occupied and so in Eqs. (2,3) we have to use the antiferromagnetic state 1$^°) of Eq.(32). In small 
clusters like the 4x4 one the one-body states are widely separated and the intra-shell interaction is 
much more important than the inter-shell one; therefore, we consider only the m states made removing 
two holes in Shf from 1$^°), neglecting the high-lying unoccupied orbitals [17]. 

We now come to item d) of Section I B, the identification of the two states which are obtained by 
spin-flip from each other and are coupled by the interaction in second order. The explicit form of 

the symmetry adapted pairs of Section V shows that {ipi^ 2 } and \4>i^ 3 } involve only k\ and k 2 . This 
suggests that by taking linear combinations we can obtain a two-body state and its spin-flipped image. 
Specifically, using the electron creation operators 



o4 ±] = ^(Ck.a ± C k2a ) 



(44) 



we build the following two-electron W = determinants: 

K) =a( +1 aJ H |0) e , \w 2 ) =a\- ] a [ + ] \0) e . (45) 

where the state |0) e is the electron-vacuum state: c^ cr |0) e = 0, Vfc,o\ With Eq.(45) we can build the 
following m states: 

\m 1 )=a [ + ] a{- ] \$ s J T F ) , \m 2 ) = a\- ] a{ +] \<P s A = F °). (46) 

Since |mi) and \m 2 ) have projection only on irreps S 2 and S3, the states (46) can be mixed by 
the operator W e ff only between themselves, and so they are the only states involved in the canon- 
ical transformation. The eigenvalue equation (12) reduces to a 2 x 2 problem; the W e ff matrix to 
diagonalize is: 



(mi|W e //|m 2 ) 
(m 2 \W e ff\mi) 



(47) 



The eigenvalues are ±(mi\W e ff\m 2 ) and the eigenvectors are ^=(1,±1). Expanding these eigen- 
vectors in the base (46), we find: 

|*1 E '> = -UK) + K» = \^ 2 )e |<&^ }, |< Ss ) = -iflmi) - |m 2 )) = |^ S3 ) e ® |*25°>, 
V2 V2 ( 48 ) 
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where the notation \ip^) e ® |$af°) stands for creating the electron-pair \ip^) over the antiferromagnetic 
state l*!*^ ) (filled shells are understood) to get a 14 hole state. We use this example to stress that 
pairing by the present mechanism is possible only if the Optimal Group Q possesses twin singlet and 
triplet representations with the same number of components, built by the same orbitals. 

In a similar way, |^> 6 ni ) an d IV^ 2 ) involve only fc 4 and fc 6 ; this suggests introducing the electron 
creation operators 

/?[±] = -L( Cfe4(7 ± Cfe6CT ); (49) 

the two-electron W = determinants 

| W ' 1 )=/3j +1 /3|- 1 |0) e , \w / 2 )=f3\- ] p [ + ] \0) e (50) 

are obtained from each other by spin-flip. Incidentally, the latter property is preserved if one switches 
to the local picture by the transformation (24): introducing a new orbital 

Citj C^u CQ a / r -i\ 

C ^ = ^T + V = 2 + V = 6 (51) 

one finds 

K> = CneTC6l|0)e, (52) 

while 

\w' 2 ) = c 6T cn a |0) e . (53) 

Thus, we have identified the m states 

\m' 1 )=p\ +] (3 [ - ] \^ s A f) , \ m ' 2 )=f3\- ] P [ + ] \<P s A =F)- (54) 

Since \m'^) and \m' 2 ) have projection only on irrcps Oi and f^, the states (54) can be mixed by the 
operator W e ff only between themselves, and we find another 2x2 matrix 



o Kiwv/K) 

(m' 2 |W e// K) 
with eigenvalues ±(m' 1 |W e //|m2) and eigenvectors ^(1,±1). Therefore, 



(55) 



jS=0\ 
- AF I 



\^) = -i=(K) + K)) = \^ ni )e ® |*I?°>, |*? 2 ) = -UK) - K)) EE |^ Q2 )e ® |^ 

V2 V2 ( 56) 

The numerical values of the eigenvalues for U — — t = 1 eV are negative in both cases, 

( mi \W eff \m 2 ) = -m.lmeV {m'^Wef f \m' 2 ) = -61.9 meV (57) 

which means that in both the cases the Cooper-like equation (12) gives singlet ground states; the 
triplets have an opposite correction to the energy. The binding energy is larger for the \ip Ql ) singlet, 
which is the true ground state (the result cannot depend on the component of the irrep). Moreover 
since l*! 5 ^ ) is a totalsymmetric singlet with vanishing total momentum, we can also predict momen- 
tum, spin and symmetry of the ground state. Our approach has enough predictive power to yield the 
symmetry, wavevector and degeneracy of the ground state. 

VIII. CONCLUSIONS 

Recently, we have proposed a new approach [15], [16] to the question of the existence and nature of 
Cooper-like pairing from repulsive interactions; this is based on a configuration mixing realized by a 
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canonical transformation; the vacuum is the ground state of the system at U = + and the Cooper-like 
bound states arise from W = pairs. Thus our approach is actually not a perturbation theory but 
its application is easier in weak coupling situations, when one can limit the configuration mixing to a 
few particle-hole excitations. The pairing mechanism involved is a form of spin-flip exchange diagram 
which is enhanced by the Ci v symmetry and our canonical transformation works independenly of the 
size of the system. We have tested the results with exact diagonalization data for open-boundary 
condition clusters [3] [17] in the three-band Hubbard model. We obtained good agreement with the 
numerical data already by the simplest approximation within our scheme, i.e. by truncating the 
configuration mixing at the level of electron-hole pairs exchange. 

In this paper we test our W — pairing mechanism within the one-band Hubbard model with 
periodic boundary condition using exact numerical data on the 4x4 square lattice [2] . In particular, 
several workers have recognised that those data could be qualitatively understood by a weak coupling 
analysis [11] [12]; however there was no clear-cut conclusion about the existence of pairing and the 
symmetry analysis of the system was not complete. 

We point out the criteria that allow one to unambiguously diagnose pairing of two holes added to 
the system. At the heart of the effective interaction there is the diagram of Figure 1 . To characterize 
the symmetry of the ground state with 2 added holes we obtained the analytic ground state at half 
filling for U = + by a new approach [21]. Moreover, we need the symmetry of the W — pairs 
in the system at hand. To this end we prove a general result that we call W — Theorem. The 
full exploitation of this theorem requires the knowledge of the Optimal Group, that is, a symmetry 
group which is so big that no accidental degeneracies arise. We have obtained this Group here for the 
first time. An adequate knowledge of the symmetry of the system allowed us to develop in detail the 
canonical trasformation on the 4x4 cluster analytically; thus we obtained good agreement with the 
data and a clear diagnosis of pairing. This result lends further support to the general approach [15], 
[16] which predicts pairing in the N x N cluster for any N; moreover, we stress that we have already 
shown how the the W = pairing grants the superconducting flux quantization [17]. 
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APPENDIX A: CONTRIBUTIONS TO THE W MATRIX FROM FILLED SHELLS 

The A 2 — body determinantal wave functions with S z = that one can build using the orbitals with 

/ 2N - 2 \ 2 

e(k) < and half of those with e(fc) = are a set of [ ^ ^ I elements. Each represents one of the 

degenerate unperturbed (U — 0) ground state configurations at half filling. First-order perturbation 
theory requires the diagonalization of the W matrix over such a basis. 

The diagonal elements of the W matrix are just expectation values over determinants \k a j kp J. . . .). 
Such an expectation value is a sum over all the possible pairs of the bielectronic elements of W like 

W(a(3, af3) = £ U{k a \n r \k a ){k \n r \k fj ) = £ U l- e ^ a -k a )r l_ e i { k,-k,)r = U_ N 2 = U_ 

the result is independent of k a and hp. Since in any determinant of the set A 2 /2 plane wave states 
are occupied for each spin, there are A 4 /4 pairs, and the diagonal elements are all equal to UN 2 /A. 
Thus, the diagonal elements shift all the eigenvalues by this fixed amount. 

The off-diagonal elements of the W matrix between determinants that differ by three or more spin- 
orbitals vanish because IF is a two-body operator. The off-diagonal elements between determinants 
that differ by one spin-orbital are sum of contributions like W(af3,jf3) — ^ r U{k a \n r \k 1 ){kp\n r \kp) 
that vanish because of the orthogonality of the plane-wave orbitals. One is left with the off-diagonal 
elements between determinants that differ by two spin-orbitals, which coincide with the corresponding 
bielectronic elements W(a(3,^S) = ^2 r U(k a \n r \ky){kp\n r \kg). This is just the matrix of W over the 
truncated Hilbert space H spanned by the states of the holes in the half filled shell, ignoring the filled 
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ones. We stress that there are N — 1 holes of each spin in Shf, thus H. is much smaller than the full 
Hilbert space of the Hubbard Hamiltonian; however, since the number of holes grows linearly with N, 
the problem is still far from trivial. 
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